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Abstract 

We show that the molecule-like configuration of DD* enables weak binding with two re- 
alistic potential models (Bhaduri and Grenoble AL1). Three-body forces may increase the 
binding and strengthen the cc diquark configuration. As a signature we propose the branching 
ratio between radiative and pionic decay. 

1 Introduction 

The motivation to study tetraquarks (also called dimesons) comes from our curiosity whether or 
not we can extrapolate our understanding of mesons and baryons (in terms of quark models) to 
two-hadron systems. Dimesons are simpler than dibaryons (or nuclear forces) since they represent 
a four-body rather than a six-body system. Heavy dimesons are cleaner than light ones since 
nonrelativistic parameterization and treatment are more justified and since they are likely to be 
longlived. Therefore in this paper we study double heavy tetraquarks as prototypes. In our extrap- 
olation from mesons and baryons to tetraquarks we assume the " A • A" type of colour dependence 
of the interaction. It seems to work properly when going from mesons to baryons. It is a challenge 
to use tetraquarks as a test whether or not this assumption is valid also for larger systems. 

However, the question arises whether or not there is a signature of tetraquarks (dimesons) 
which could clearly distinguish their decay from the decay of two independent mesons. Consider 
the analogy. The mass of the neutron is 1.3 MeV larger than the mass of the proton, making the 
neutron unstable against the weak interaction and results in the n — > pe~i7 e decay. But when the 
neutron is bound in the deuteron with a binding energy of -2.23 MeV, the decay is kincmatically 
forbidden and the neutron becomes stable. If we replace the two baryons in the deuteron with two 
mesons we obtain the dimeson. When one of the mesons is a vector meson, its dominant decay 
mode in the case of weak binding would be the radiative decay B* — > B^/, or in the system of 
the D mesons D* — > Dj, as well as the strong decay D* — > Dir. The B*B dimeson, however, is 
probably bound strongly enough so that radiative decay becomes energetically forbidden [1, 3, 5] 
and it can decay only weakly. The binding of the D*D dimeson is expected to be much weaker, 
but it might still be stable against strong interaction so as to decay only electromagnetically, or, 
at least, the strong decay might be considerably suppressed. The dramatic change of the decay 
modes and lifetime of vector mesons can serve as a tool for detecting tetraquarks with nonzero total 
spin. This effect would be very helpful in situations where the binding energy is small compared 
to experimental errors so that the detection of the tetraquark from the invariant mass of final 
particles could not distinguish between events where the two initial mesons were either free or 
weakly bound before the decay. 

In this study we consider tetraquarks made up of the same type of mesons, namely DD* and 
BB* tetraquarks with quantum numbers I — 0, S — 1, P — +1, since they are the best candidates 
for binding with respect to the DD* (BB*) threshold, as was already noted previously [1, 2, 3, 4, 5]. 
Such states cannot decay strongly or electromagnetically into two B or two D mesons in the S 
wave due to angular momentum conservation nor in the P wave due to parity conservation. 



There are two extreme spatial configurations of quarks in a tetraquark. The first configuration 
which we call atomic is similar to A;,, with a compact bb diquark instead of b, around which the light 
antiquarks are in a similar state as in the Ah baryon. The stability of the double heavy tetraquarks 
with this structure was first investigated by Lipkin [6]. The second configuration which we call 
molecular is deuteron-like, where two heavy quarks are well separated and the light antiquarks are 
bound to them similarly as in the case of free mesons. This configuration is more likely to appear 
in weakly bound systems and was studied by Manohar and Wise [7] and Tornqvist [8, 9] in the 
framework of the pion exchange between two heavy mesons. It has been shown that nonrclativistic 
potential models in general give rise to atomic structure for the BB* tetraquark [1, 3], due to the 
large quark mass asymmetry. We use the BB* system as a benchmark against which we compare 
other tetraquarks. 

In Sect. 2 we repeat calculations of the BB* tetraquark with the Bhaduri potential [10] and 
also with the AL1 potential [11] which due to the additional mass-dependent smearing of the spin- 
spin interaction gives a better description of meson spectroscopy. We then present results for the 
DD* tetraquark, which is on the verge of being either bound or a resonant state, depending on the 
effective potential used in the calculations. We show that its structure is molecular, but with the 
introduction of the three-body force (Sect. 3) it can become atomic. The fact that this system is 
so close to the D + D* threshold makes it very sensitive to the details of the effective interaction 
and therefore a promising candidate for studying the nature of the effective interaction between 
constituent quarks in nonrclativistic potential models. The estimated production rate is not high 
but tolerable. Due to the strong influence of weak binding on the decay channels, it presents a 
very interesting experimental situation (Sect. 4). 

2 Bound states of heavy tetraquarks 

The general idea of possible stable heavy tetraquarks has been first suggested by Jaffe [12]. It soon 
became clear that the systems with unequal masses of the quarks and the antiquarks are more 
promising, since the binding energy strongly depend on the mass ratio [13, 14, 15]. The b/u mass 
ratio is shown to be large enough to make Tbb stable, while the c/u mass ratio is under-critical 
for atomic structure. The calculation of various tetraquarks in the harmonic oscillator basis [3, 4] 
have shown that only two tetraquark systems have their energy lower than the two-free-meson 
threshold, namely bbud (1=0, J=l) which we denote as Ty, and bbsu or bbsd (1=1/2, J=l), while 
besu or bbsd (1=1/2, J=l) lie on the threshold. For deeply bound states, these results should be 
very accurate but since this basis cannot accommodate asymptotic states of two free mesons, there 
is an open question whether or not weakly bound states of two mesons have been missed. 

In our work we use an expansion in the basis proposed by [5] which is described in Appendix 
A. There are three different sets of internal coordinates. The first one (Fig. 10 a)) is convenient for 
expansion of those strongly correlated and deeply bound tetraquarks where we expect the atomic 
structure in which the diquark in Tqq formed by two heavy quarks plays a similar role as the heavy 
b antiquark in the Ah baryon, while the light quarks in both systems are in the same radial, spin, 
colour and isospin configurations. The second and third sets from Fig. 10 represent the direct and 
exchange meson-meson channels. These configurations are needed to build up the basis for the 
two free mesons - the threshold state, and are also of crucial importance for searching for weakly 
bound tetraquarks where the molecular structure would be dominant. 

We search for eigenstates of our Hamiltonian using the variational method, applying a general 
diagonalization of the Hamiltonian (Appendix B) spanned by the non-orthogonal basis functions 
constructed in Appendix A. We built the basis functions step by step by adding the best config- 
urations from Fig. 10 with the best colour-spin configurations allowed for our quantum numbers 
(IS=01, positive parity and colour singlet), after having optimized the corresponding Gaussian 
widths. In order to obtain a 0.1 MeV accuracy we constructed bases in this way with up to 
N m ax = 90 and N max = 140 functions for the Thh and T cc tetraquarks, respectively (Appendix D). 
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These basis states can also accommodate two asymptotically free mesons if the four-body problem 
has no bound state. 

2.1 T bb 

First we test our method on the Tbb system and compare our results with Ref. [3, 4]. In our 
calculations we use two one-gluon exchange potentials, the Bhaduri and AL1 potential. Their 
properties are described in Appendix B. The Bhaduri potential quite successfully describes the 
spectroscopy of the meson, as well as baryon ground states. This is an important condition since 
in the tetraquarks we have both quark-quark and quark- antiquark interactions. The AL1 potential 
is just an improvement of the Bhaduri potential where the smearing of the colour-magnetic term 
in the Hamiltonian depends on the masses of the quarks. This then results in better quality of the 
meson spectra, in particular in the charmonium sector, where the Bhaduri potential predict only 
half of the observed hyperfinc splitting between the rj c and J/ip state. 

The results obtained with the Bhaduri potential are presented in the first three columns of 
Table 1 where they are compared with results from [3]. In the last three columns our and [4] 
results for the AL1 potential are listed. Since the harmonic oscillator basis cannot accommodate 
two asymptotically free mesons, one obtains a positive binding energy, as for example for spin 1, 
isospin 1 state calculated with the Bhaduri potential. 

Table 1: The mass of the Tbb tetraquark. Column 1: spin S, isospin I, Column 2: lowest meson- 
meson threshold for the Bhaduri potential, Column 3: our results, Column 4: results of ref. [3] where 
expansion in the harmonic oscillator basis was used, Column 5: lowest meson-meson threshold for 
the AL1 potential, Column 6: our results with the AL1 potential, Column 7: results of ref. [4] 
where the same basis was used as in ref. [3]. 
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Figure 1: Probability density of two heavy quarks pbb, of two light antiquarks p m and of a light 
antiquark and heavy quark pqb in Tbb as a function of interquark distances for the AL1 potential 
(q is ii or d). 
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To obtain a better understanding of the Tbb tetraquark we now turn to the radial, spin and 
colour structure of this system. The probability densities for finding quark i and (anti)quark j at 
the interquark distance as shown in Fig. 1 is calculated via 

Pij(r) = (ip\S(r - rij)\i>). 
The projection of this probability density on the colour triplet states |3i2334)c is 

pf^-\r) = W>|3 12 3 3 4}c(3i 2 3 3 4|c^ - r y )|V>>, 

and similarly for the other spin and colour projections presented in Fig. 2. 

From Fig. 1 it can be seen that the probability density of two heavy b quarks is strongly 
localized so one can make a rough approximation in which the heavy diquark is pointlike. The 
average distances between a light antiquark and a heavy quark, and between a light antiquark and 
a light antiquark are almost the same, so the heavy diquark and the two light antiquarks form 
some sort of equilateral triangle which can also be seen as a "Y" shape configuration. 

In Fig. 2a we see that the dominant colour configuration is 3i2334 where two heavy quarks are 
in the colour antitriplet state. The 612634 configuration becomes relatively important at very large 
distances where the absolute probability density is negligible. The ratio of these two configurations 
approaches 2 for large separation between two b quarks which is consistent with the fact that for 
large distances we have two white mesons and the octet configuration is not present. This can 
be seen from the decomposition of colour octet and singlet (quar k- antiquark)- (quark-antiquark) 
states into colour sextet and triplet diquark-antidiquark states. 

|li3l24)c 
|8i38 2 4)c 

|ll4l 2 3)c 
(814823)0 

Since the heavy diquark is in a spatial symmetric and colour antitriplet state, it must be due to 
the Pauli principle in the spin symmetric S=l state. This can also be seen in Fig. 2b. In this 
figure the |li2l3 4 )s configurations are not shown since they are three orders of magnitude smaller 
and can thus be neglected. 

We have shown that the Tbb tetraquark has an atomic Aft-like structure where the heavy diquark 
in the colour antitriplet state can be approximated with a heavy antiquark b, while the light anti- 
quarks are in isospin and spin state just like in Aft baryons. This then justifies the assumptions 
made in [1] that the light quarks in Tbb have a similar spatial function as in Aft. The binding energy 
can then be phenomenologically estimated [1] to be -134 MeV for the AL1 potential and -139 McV 
for the Bhaduri potential, which is close to the detailed calculations presented here. This result is 
independent of the form of the interaction between u and d since if the ud wave function in T fcb 
and Aft is the same, then also the interaction matrix element is the same in both cases and cancels 
out in the comparison. For this reason, Xftft with its Aft-like structure cannot distinguish between 
one-gluon (OGE) [10, 11] and Goldstone boson exchange (GBE) [16]. As argued in Appendix C, 
it is crucial to have correct Aft mass to ensure the correct ud interaction and the GBE calculations 
which do not fit the Aft mass [17, 18] may strongly overbind the Tbb as well as T cc . 

Because of a dominant colour triplet-triplet and spacial "Y" shape structure of the tetraquark, 
the introduction of a weak colour dependent three-body interaction would only shift the mass 
of the tetraquark and not produce any significant changes in the wave function, similarly as in 
the baryon sector. This effect can be compensated by a reparametrisation of parameters in the 
two-body potential or constituent quark masses so as to reproduce correct baryon ground states; 



= y^|3i 2 3 3 4)c + y||6i26 3 4)c, 
= -^/||3i 2 334)c + ^/||6i 2 6 3 4)c- 
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Figure 2: Results of calculations with the AL1 potential, a): Probability densities p^ lp '^ and 
Pbb* 1 ^ OI tne two neavv quarks projected on the colour triplet |3i2334)c and colour sextet states 
|3i2334)c, respectively, and their ratio, b): Probability density and p^P? of the two heavy 
quarks projected on the spin 1 states |li2034)s and spin states |0i2l34)s, respectively. 



the shrinking of baryon spectra due to three-body interaction does not concern our considerations. 
Therefore the T^b tetraquark is unsuitable for studying the influence of the three-body interaction. 

There are also experimental problems with the T^b tetraquark. The only promising production 
mechanism - double two-gluon fusion (g + g — > b + b) 2 gives a very small production rate of about 
5 events per hour at LHC [19], [20]. Moreover, since T&j is below the total mass of two B mesons, 
it can decay only weakly and thus has no characteristic decay different from separate B decays. 
This then also makes the T^b tetraquark unpromising from the experimental point of view. 

2.2 T cc 

The T cc tetraquark is much more promising than the Tbb tetraquark. It can be more easily produced 
and detected (Sect. 4) and we shall see that it better discriminates between different binding 
mechanisms. 

With the expansion of the tetraquark wave function in the harmonic oscillator basis one cannot 
find any bound state for the T cc system with the Bhaduri or AL1 potential. But as mentioned in 
the previous subsection, this can also be due to the fact that this method can miss weakly bound 
states. And this is exactly what happened as one can see from Table 2 where our results are 
presented. With both potentials, a weakly bound state does appear. 



Table 2: The mass of the T cc (S=l, 1=0) tetraquark. Column 1: type of potential, Column 2: 
lowest meson-meson threshold for a given potential, Column 3: our results, Column 4: results of 
ref.[3, 4] where expansion in the harmonic oscillator basis was used. 
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Figure 3: Probability density of the two heavy quarks p cc , of the two light antiquarks pqq and of a 
light antiquark and a heavy quark pq c in T cc as a function of the interquark distance, a): results 
for the AL1 potential, b): results for the Bhaduri potential 



As in the previous subsection for the Tbb we now repeat the two-quark probability density 
analysis for the T cc system. In Figs. 3 and 4 the probability densities their projections on colour 
and spin states as a function of interquark distance are shown. 




r (fm) r (fm) 



Figure 4: Results of calculations with the AL1 potential, a): Probability densities pL" and 
Pcc eXt ^ of the two heavy quarks as a function of interquark distance projected on the colour triplet 
state |3i2334)c and on the colour sextet state |3i2334)c, respectively, b): Probability densities 
pec^ ) Pcc^ and pi} for the two heavy quarks as a function of interquark distance projected on the 
spin 1 states |li2034)s and |li2 134)5 and on spin state |0i2 134)5, respectively. 

We can see in Fig. 3 that the wave function between heavy quarks is much broader and has 
an exponential tail at large distances. If we look at the structure of the quark-quark probability 
density in Fig. 4a we see that at around r <~ 1 fm sextet configurations become larger than triplet 
ones and soon after the ratio of the colour configurations stabilizes at 2. This supports the picture 
of molecular binding of the D and D* meson in the T cc . This can also be confirmed from the 
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results shown in Fig. 4b where at distances larger than 1 fm the probability for two heavy quarks 
in spin and spin 1 states is equal. This follows from spin recoupling 



|ll3024)s = ~\[\ |ll2l34)s + ^75(|ll2034)s - |0l2l34)5) 
|0l3l24),S = ~\J\ |ll2l34)s — ^(|ll2034)s — |0i2l34)s) 



If we assume that in the ground state orbital and also colour wave function is symmetric/ antisymmetric 
with respect to permutation of identical particles (~ |li3l-24)c + 1 ~ |li4l23)c)i the spin wave func- 
tion must be antisymmetric and has the form 

|V) ~ l/\/2(|li 3 24 )s - |0i 3 l 24 )s) - -1/V2(|li 2 34 )s - |0i 2 l 34 )s). 

The contribution of the |li20 3 4)s and of the |0i2l-34)s configurations is thus equal. Similar con- 
clusion also holds, if we recouple to |li4023)s and |0i4l23)s- All relevant spin recouplings can be 
found in [5]. 



PccW 






Figure 5: Decomposition of the quark-antiquark probability density pq c (bottom right) into contri- 
butions corresponding to different distances between heavy quarks (left). The details are explained 
in the text. 

The probability density of the light antiquarks pqq shows similar behaviour, while the radial 
dependence of the quark-antiquark probability density pq c has a strong peak at smaller distances 
and a twice lower tail than the pqq and p cc (Fig. 3). How this structure appears is shown in Fig. 5. 
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Due to symmctrization of the configurations into which we expand the tetraquark wave function, 
the first light antiquark is bound to both heavy quarks ci and C2 with equal probability. The 
probability density Pq C (r,r cc ) of finding the first light antiquark at the interquark distance r from 
the first heavy quark when the distance between the heavy quarks r cc is r\ , and is shown on 
the left hand side of Fig. 5. The heights pi are proportional to the probability density of finding 
two heavy quarks at the interquark distance r i; as is schematically depicted on the top right of 
Fig. 5. After integrating over all possible interquark positions j-j we obtain a strong peak when 
the light antiquark is bound to the first heavy quark c\ and a long tail when it is bound to the 
second heavy quark c^. If we ignore the interference between these two situations which appears 
when ri is smaller than the size of the free D meson, half of the probability is in the peak and the 
other half in the tail. 

The colour and spin decomposition of the quark-antiquark probability density pq c is shown in 
Fig. 6. We see in Fig. 6b that at small distances (r < 1 fm) where the peak is situated, the 
probability density between first heavy quark (particle 1 from Fig. 10) and first light antiquark 
(particle 3 from Fig. 10) corresponds to the pseudoscalar D meson |0i3l24)s and the vector D* 
meson |li3024)s- To understand the spin structure of the tail, we also present in the same graph 
the projection on the spin |li4023)s, |0i4l23)s an d |li4l23)s states. We see that the dominant 
spin |li3l24)s channel here can be decomposed into the |0i4l23)s and |li4023)s contribution. This 
means that for larger separations r the heavy quark dresses with the second light antiquark (particle 
4) into the D or D* meson. From this we can conclude that there is no appreciable contribution 
from the D*D* configuration also for large r. For large r on can see in Fig. 6a that the first 
light antiquark combines with the second heavy quark into a colour singlet state. The octet colour 
dominance at small interquark distances is due to the formation of diquark-antidiquark structure, 
similar as to that in Tbb- But in T cc this is not a dominant structure, it represents only about 
a third of the total probability in the case of the AL1 potential and even less for the Bhaduri 
potential. 




r [fm] r [fm] 



Figure 6: Results of calculations with the Bhaduri potential, a): Probability densities pf^ 
and p^° ct '' 1 of the quark and antiquark projected on the colour singlet |li4l23)c and colour octet 
states 1814823)07 respectively, b): Projection of Probability density pq c on various spin states. 
The |0i3l24)s and |li3024)s projections are almost exactly (after renormalization) the probability 
densities of D and D* mesons, respectively. 
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We also analysed the decrease of the tail of the probability densities from Fig. 3a. At large 
distances one can approximate the tetraquark as a bound state of two mesons which is described 
by the Schodinger equation 



2m r dr 2 



t/j + V(r)ip + Ei/>, 



where m r is the reduced mass of D and D*. Assuming that at large distances the colour wave 
function is singlet, the potential V(r) should approach zero and thus the asymptotic behaviour of 
V> should be 



ip(r — > 0) — > exp(-Kr), k = \/\E hind .\M Icd Jhc, 



(1) 



where -Ebind. is the binding energy of the tetraquark (2.7 MeV for AL1) and M re( j. the reduced mass 
of the D and D* meson. This can be very clearly seen in Fig. 7 where we plotted the logarithm 
of the wave functions. The quark-antiquark probability density is multiplied by 2 as explained in 
the previous paragraph. All three interquark probability densities follow the predicted exponential 
decrease from Eq. 1. 
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Figure 7: Logarithms of probability density as a function of interquark distances compared with 
the analytically expected slope 



3 Three-body interaction 

The idea of introducing a small amount of three-body interaction into nonrelativistic potential 
models is not new. In the baryon sector this additional interaction can be used to better reproduce 
the ground state spectroscopy [11]. But generally a slight three-body interaction cannot produce 
any significant changes in other baryon properties, and the desirable shift of levels can also be 
reproduced by modification of the parameters in the two-body potential. 

Before we introduce such an interaction into the tetraquark system we shortly discuss the 
structure of this interaction. For the radial part we take the simplest possible radial dependence - 
the smeared delta function of the coordinates of the three interacting particles. Since we are working 
with the colour-colour type potential it is natural that also the three-body potential possesses some 
colour structure. The colour factor in the two-body Bhaduri or AL1 potential is proportional to the 
first (quadratic) Casimir operator ; C^ 1 ^ = A • A. Following this, we introduce in the three-body 
potential the second (cubic) Casimir operator = d abc X a ■ A& • A c . A deeper discussion of 
the properties that the colour dependent three-body interaction must fulfil can be found in [26], 
[27], [28]. 
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It should be noted that in the baryon sector such a colour structure is irrelevant since there is 
only one colour singlet state and thus the colour factor is just a constant which can be included into 
the strength of the potential. In tetraquarks the situation is different since there are two colour 
singlet states: §12334 and 612634 (or I13I24 and 813824 after recoupling). The three-body force 
operates differently on these two states and one can anticipate that in the case of the weak binding 
it can produce large changes in the structure of the tetraquark. This cannot be otherwise produced 
simply by reparameterization of the two-body potential, so the weakly bound tetraquarks are a 
very important laboratory for studying the effect of such an interaction. 

The form of the three-body interaction we introduced into the tetraquark is 

V£ 9 {fi,?j,?k) = -Id^X-X^UocM-irl+r^+rl)/^}, 
V q %(n,fj,fk) = \d abc \«\ b *\l*U eM-(rl+r%+rl)/rl\. 



Here m is the distance between i-th and j-th (anti)quark, and similarly for rjk and rfej. A a are 
Mann colour matrices and d abc arc the SU(3) structure constants ({A a , A 6 } = 2d abc X c ). 
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Figure 8: a): Probability density between two c quarks in the T cc tetraquark as a function of 
interquark distance for three different values of the strength of the three-body potential, b): Mass 
of the T cc tetraquark as a function of the smearing parameter ro in the three-body interaction for 
four different strengths Uq of the three-body potential. 

The diagonal matrix elements of the colour part of the three-body interaction between two 
quarks and an antiquark are -5/18 and 5/9 for IS12334) and I612634) colour states, respectively. 
If the strength of this interaction Uo is negative it will lower the states with diquark-antidiquark 
configuration and increase the binding as can be seen in Fig. 8b. For a strong three-body interaction 
T cc loses the molecular structure and the triplet-triplet colour configurations become dominant and 
the T cc tetraquark becomes similar to Tbb- A drastic change in the width of the probability density 
can already be seen for strength Uo = -20 MeV. In the baryon sector such an interaction would 
merely lower the states by about Uo so it would have no dramatic effect nor would it spoil the fit to 
experimental data. Since the predicted energies of ground state baryons for the Bhaduri and AL1 
potential are above the experimental values, this is actually a desirable feature. The dependence 
of the mass of the T cc tetraquark on the strength of the potential Uo and on the smearing of this 
potential is shown in Fig. 8b. One can see that for large smearing the mass of the tetraquark is 
shifted by about 4yg[/o in agreement with the dominance of the triplet-triplet colour configuration 
in this state. 
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4 Production and detection of T, 



We shall focus only on the possibility of detecting the T cc tetraquark, since even at LHC the 
production of the tetraquark would be below the rate where one could hope to detect it. 

Since double charmed baryons were probably detected at SELEX [21] one can expect that if 
the T cc is bound, it was also produced there, with a production rate about ten times smaller than 
double charmed baryons. This estimate is based on the fact that a heavy quark gets dressed with 
a light antiquark into a heavy meson with a probability of roughly 0.9, and with a probability 
of 0.1 it combines with two light quarks into a A baryon [22]. But since SELEX found, with 
their cuts, only fifty candidates for double charmed baryons, the statistics for detecting double 
charmed tctraquarks should be improved. Another experiment where one might look for T cc could 
be provided by LHC where one can estimate the production rate as large as 10 4 events/hour [23]. 

However, the T cc tetraquark has a molecular structure in which the mesonic wave function is 
not strongly influenced by the other meson. The large mean square radius of the T cc tetraquark 
has consequences for the production mechanism. It is no longer required that the two c quarks 
come close to produce first a diquark which then gets dressed by light antiquarks. This can give a 
larger production rate at SELEX than estimated above. This also makes machines like the RHIC 
ion collider [24] interesting candidates for searching for tctraquarks. 

There is also an interesting possibility of production and detection of the T cc tetraquark in B- 
factories. Belle [29] has reported a measurement of double charm production in e + e~ annihilation 
at y/s — 10.6 GeV and found that 

cr(e + e~ — > J/ipcc) <~ lpb, 

which corresponds to about 2000 events. Since the total mass of four D mesons is close to the cm. 
energy, the c quarks created in this process have small relative momentum which is very important 
in the T cc production. 

The main problem with detection of the weakly bound T cc tetraquark is how to distinguish 
the pion or photon emitted by the decay of the free D* meson from the one emitted by the D* 
meson bound inside the tetraquark. We can exploit the fact that the phase space for D* — > D + tt 
decay is very small. Note that m D »+ — m D + — m^o — 5.6 ± 0.1 MeV, mrj*o — m D o — m n o = 
7.1 ±0.1 MeV, m D »+ — to d o — — 5.87± 0.02 MeV. This has a strong impact on the branching 
ratio between radiative and hadronic decay. Since the D* meson inside the tetraquark is not 
significantly influenced by the other D meson (Fig. 5) in the tetraquark, we expect that the partial 
width for the magnetic dipole Ml transition would be very close to the width of the free meson 
while the width for hadronic D* — > D + ir decay will decrease with stronger binding and will become 
energetically forbidden below the D + tt threshold. The hadronic decay of the T cc tetraquark is a 
three-body decay which is commonly represented by the Dalitz plot. 

Let us assume that the T cc tetraquark is below the D + D* threshold but above the D + D + 7 
and D + D + ir as was the case in our nonrelativistic potential models. Then the partial decay rate 
for the T cc — >D+D+7r is given by 

dT = ^^JWdm^dmi, (2) 

where particles 1 and 2 are two final D mesons and particle 3 is a tt emerging from the decaying 
tetraquark. Here m\ 2 — (pd + Pd) 2 and m?, 3 = (pd + Pn) 2 and M is the mass of the tetraquark. 
Since the total masses of the D* + D and 2D + tt are so close there is a strong isospin violation in 
the decay which cannot be reproduced with the Bhaduri or AL1 potential where the D* and the 
D isospin doublets are degenerate. We shall not try to modify the interaction to accommodate 
the dependence of the decay on the isospin of the particles, but we shall rather work with the 
experimental masses taken from the PDG [30]. The allowed region of integration over dm\ 2 and 
dm 2 3 for three different binding energies is plotted in Fig. 9a. If we assume |.M| 2 is constant, which 
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is very plausible in our case, the allowed region will be uniformly populated with experimental 
events so that the measured partial decay rate T will be proportional to the kinematically allowed 
area from Fig. 9a. The dependence of this area I as a function of the binding energy of the 
tetraquark is shown in Fig. 9b. 

Up to now, we have discussed only true bound states, but the T cc tetraquark can also be a 
resonant state above the D + D* threshold. Then if the resonance is situated near the threshold, 
beside the T cc — > D + D* decay there will still be a significant fraction of hadronic T cc — > D + D + n 
decays. This region of positive binding energy is not presented in Fig. 9. But one can see that in 
a similar manner as in the case of weakly bound tetraquarks the low-lying resonant state can be 
identified from the Dalitz plot. 

Table 3: Mean distance between two heavy quarks (r cc ) in the T cc tetraquark and the value of the 
centrifugal potential for L=2 state for the Bhaduri and the AL1 potential. 





( r cc) 


(ti z L{L + l)/2M red )(l/rl c ) 


Bhaduri (L=2) 


2.4 fm 


174 MeV 


AL1 (L=2) 


1.6 fm 


232 MeV 



As a remark, we present an alternative way of detecting the weakly bound tetraquark. Since 
the mean radius of the tetraquark is large, as one can see from Table 3, the centrifugal barrier 
for the L=0 — > L=2 transition is comparable with the available energy in the D* decay, so there 
is also a possibility of the electric quadrupole transition E2. This is a two pion exchange process 
which is beyond the scope of the potential model used here. 
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Figure 9: a): Dalitz plot, b): Integrated Dalitz plot. 



5 Conclusions 

We have shown that in popular nonrelativistic potential models the T cc tetraquark is bound against 
the DD* threshold and that it has a molecular structure. Therefore the approximation based on 
the assumption of the atomic Afc-like structure which suggests that the system is not bound [1] is 
not valid. This dramatically different situation as compared with the Tbb tetraquark makes the T cc 
tetraquark an interesting laboratory for more profound studies of the nature of the interactions 
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between quarks, since such a weakly bound system is more sensitive to the detailed features of the 
interaction. 

As a signature for the T cc — DD* tctraquark one might exploit the very small phase space of 
the D* — > Dtt decay which is very sensitive to the binding energy of D* to D. 

It seems tempting to compare the DD* dimeson with the recently found DD* (3.872) state 
which is just above the DD* threshold. This is, However, a completely different situation due to 
the low J/ipir or J/iprj thresholds, which arc degenerate for both potentials used here. Thus a 
delicate coupled channel calculation would be needed in the search for resonances which is beyond 
the scope of this paper. Our energy minimisation procedure gave in fact the threshold energy 3232 
McV for spin 1 ecu d state calculated with the Bhaduri potential, consistent with ref.[3], where due 
to their basis they obtained somewhat higher value of 3468 MeV. 

We did not study in detail the ccsu tetraquark since it has a lower threshold, D S D rather than 
D S D* or D*D, and it is not likely to be bound. 

Acknowledgement. The authors would like to acknowledge the encouragement by Jean-Marc 
Richard to look whether or not the DD* system might be bound after all. Acknowledgment is also 
due to the collaboration with Danielle Treleani and Alessio Del Fabbro on questions of production 
of the Ttb and T cc tetraquarks. 

This work was supported by the Ministry of Education, Science and Sport of the Republic of 
Slovenia. 

A Configurations 

We express the orbital part of the tetraquark wave function in terms of Gaussians. The coordinate 
systems used here are shown in Fig. 10. The transformation between various coordinate systems 
and some details about the calculation of the kinetic and potential matrix elements are given in 
[5]. 




Figure 10: Two quarks (dashed circles) and two antiquarks (empty circles) in three different relative 
coordinate systems. The orbital wave function is then a Gaussian function of relative coordinates, 
a) diquark-antidiquark: K\{Cij) — exp(— XjCyXj), b) direct channel: K 2 {Cij) = exp(— YiCijyj), 
c) exchange channel: Ks(Cij) — exp(— ZidjZj). 

The most general form for the ground state of the tetraquark (L=0) can be expanded in any 
of the three configurations given in Fig. 10 

R = J2CnK r (C?j), CI = C£, r = lor2or3. 

n 

So far it looks as if not all coordinate systems were needed, but in numerical calculation it is 
convenient to limit the test functions to those in which CV, = if i ^ j. With this we reduce our 
problem of optimization of 6 parameters Cij into three optimizations (r = 1, 2, 3) of 3 parameters 
Cu — {ci, 02,03}. Though this somewhat restricts our Hilbert space we still expect that it would 
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not have any significant effect on the calculation of the ground state, because by using all three 
configurations from Fig. 10 we can have nonzero relative angular momentum between two quarks 
?12 or two antiquarks £ 34 by using the systems b) and c) but still keep the total angular momentum 
zero. The orbital part of the wave function has thus the form 



R = Y J C n K r {c n 1 ,c n 2 ,cV i ), r=l,2,3. 



The tetraquark wave function must posses correct symmetry against the permutation of the 
two quarks P(l, 2) or antiquarks P(3, 4). The effect of these permutations on relative coordinates 
are 



P(l,2) 
P(l,2) 
P(l,2) 



Ki(ci)] =Ki(ci) 
K 2 (ci,C2, c 3 ) 

K 3 (ci,C 2 ,C 3 ) 



P(3,4) 

= K 3 {c 2 ,c 1 ,c 3 ), P(3,4) 
= K 2 {c2, Cl ,c 3 ), P(3,4) 



K 2 { Cl ) 
K 3 { Cl ) 



= K 1 (c i ), 
= K 3 (a), 
= K 2 {ci). 



We see that if c\ 7^ c 2 the configuration K 2 and K 3 loose symmetry properties with respect to 
permutations of the heavy quarks. To obtain functions with good permutation symmetry we 
make linear combinations of configurations. For shorter notation we introduce K j(c2, C\, c 3 ) = 

K i {c 1 ,C2,C S ) 



Ri(ci) =Jfi(c 
R 2 (c t ) =K 2 (c 

ifefc) =K 2 {c, 

Ri{Ci) =K2(C, 

R 5 {ci) =K 2 (c, 



). 

)+X 3 ( Cj ) + if 2 (c J ) + ^ 3 (c J ), 

)-Jf 3 (Ci)+^2(Ci)-^ 3 (Ci), 

)-K 3 (ci)-K 2 (ci)+K 3 {ci), 
) + K 3 ( Ci ) - K 2 (ci) - K 3 ( Ci ). 



The effect of permutation of identical (anti)quarks is then 

P(12)[i2i(c 
P(12)[P 2 (c, 
P(12)[P 3 (c, 
P(12)[P 4 (C 
P(12)[P 5 (c, 

In the spin space of four quarks we have three different spin 1 representations. Most suitable 
basis for studying permutation properties is obtained by coupling the quarks into a diquark and 
antiquarks into a antidiquark. The three basis states are then 



)] 


= P(34)[i2i(ci) 


= Ri(ci), 


)] 


= P{M)[R2{ci) 


= R 2 (c l ), 


)] 


= P(34)[i? 3 (c 4 ) 


=-Ra(ci), 


)] 


=-P(34)[P 4 (cj) 


= Ri(ci), 


)] 


=-P(34)[P 5 (ci) 


= -R 5 (ci) 



I ll2 , O34) , |0i2 , I34) , |ll2, I34), 

and the permutation of identical particles on this states gives 



(3) 



P(12)[|li2,034>] = |ll2,0 34 ), 
P(12)[|0i2,l 34 )] - — |0i2, I34) 
P(12)[|l 12 ,l 34 )] - |ll2,l 34 ), 



P(34)[|l 12 ,0 34 )] = -|li2,0 34 ), 
P(34)[|0 12 ,l 34 }] - |0i2,l 34 ), 

P(34)[|l 12 ,l 34 }] - |ll2,l34>. 



In the colour space there are two different colour singlet representations which in the diquark- 
antidiquark basis can be expressed as 



|3i2,3 34 ), 1612,634), 



(4) 
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and wc have 



P(12)[|3 12 ,3 34 }] = -|3 12 ,3 3 4>, P(34)[|3 12 ,3 34 )] = -|3 12 ,3 34 ), 
P(12)[|6i 2) 6 34 >] = |6i 2) 6 34 ), P(34)[|6i2, 634)] = |6 12 , 6 34 ). 



From the original three spatial configuration shown in Fig. 10, the three spin configurations 
given in Eq.(3) and two colour configuration of Eq.(4) one can build up 3 • 3 • 2 = 18 function. 
Eight of them are antisymmetric with respect to exchange of heavy quarks and symmetric (isospin 
0) with respect to exchange of light antiquarks. 

V>i = J?i(ci)|3i 2 3 34 ) c |li20 34 )s, 1P2 = -Ri(ci)|6i 2 6 34 )c|0i 2 l 34 )s, 

1p3 = #2(Ci)|3i 2 3 34 )c|ll20 3 4)s 5 ^ 4 = P2(Ci)|6i 2 6 34 )c|0i 2 l 34 )s, 
1p5 = P3(Ci)|6i 2 6 34 )c|ll20 34 )s, i>6 = P3(ci)|3l 2 3 34 ) C |0i 2 l 34 )s, 
1p 7 = i? 4 (c l )|3i 2 3 34 )c|ll 2 l 34 )5, i>8 = P5(Ci)|6l 2 6 34 )c|ll 2 l 34 )s- 



For better description of weakly bound states we also add additional configurations which cannot 
be decomposed into a simple product of orbital, colour and spin parts. 

^9 = ((K 2 {Ci) + K 2 (Ci))\l 13 l 24 ) C + (K 3 (d) + ^ 3 (Cj))|ll4l23>c) |0i 2 l 34 )s, 

VlO = (^(^2(c J )+^2(c i ))|ll 3 l 24 >C- (^3(c i ) + X 3 (c 4 ))|l 14 l 23 )c^)|ll 2 34 )5, 

V'll = (K 2 (ci) -^2(Cj))|ll 3 l24)c|0l 3 l24>S+ (K 3 (d) - ^(Ci)) |li 4 l23>c|0l4l24)s, 
= (K 2 { Ci ) - ^2(Cj))|ll 3 l24>c|ll3024>S + (K 3 {c t ) - K 3 (Ci)) | ll 4 l 23 )c | ll 4 23 )s, 
^13 - (K 2 ( Cl ) - £2(Cj))|ll3l24>c|ll3l24>S + (^s(Ci) - K 3 (<*)) | ll4l 23 >C | ll 4 l 2 3>S. 

It is obvious that these configurations also respect permutation symmetry. 

If we have a strong quark mass asymmetry we expect clustering of heavy quarks into a diquark 
which then results in the atomic structure of the tetraquark [5] , so that the first coordinate system 
in Fig. 10 is more suitable and the dominant colour configuration has the diquark in antitriplet and 
the antidiquark in triplet colour state. On the other hand, if the binding is weak, the direct and 
exchange meson-meson channels are more adequate and the important configuration has singlet- 
singlet colour structure. 



B Potential models 



For solving Schrodinger equation in non-diagonal basis we use general diagonalization of the Hamil- 
tonian 



=1. .3; 
= i + l. .4 



fe, i = 1, ..13; m,n=l,..N max ; 

where is the m-th basis function (See Appendix D) and N max is the dimension of the basis. 
The kinetic energy operator written in the basis a) of Fig. 10 has the form 

Wk\W k \W) = -6(Vm n >Tr[(C m + C n Y x C m TC n 



T = 



h 2 c 2 





mj — mi 
2v / 2mim2 





m 3 +m 4 
277137714 
m^—rri4 
\m3m4 



2V2: 



2y/2mim2 
7713—7714 

2V2wi3m4 
1 y^4 j_ 
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We tested two different potentials. First one was proposed by Bhaduri and collaborators [10] 
and the improved one proposed by Silvestre-Brac and Semay [11] which we will denoted by AL1 
potential. 

- Bhaduri potential: 

,rB Af Xff TT a a h 2 e-^/r° \ 



Tin Vi Tj 



m b = 5259 MeV, m c = 1870 MeV, 

m s = 600 MeV, m u = m d = 337 MeV, 

U = 685 MeV, a = 77 MeVfm, 

/3 = 706.95 McV/fm, r = 0.4545 fm. 

AL1 potential 



B 



r (m,i,mj) = A 



rrii + m 



01 

iiu-j 



m b = 5227 MeV, m c = 1836 MeV, A = 1.6553 GeV B_1 , 

m s = 577 MeV, m u = m d = 315 MeV, B = 0.2204, 

U = 624.075 MeV, a = 74.895 MeVfm, 

(3 = 629.315 MeV/fm, a = 274.948 MeVfm. 



C Effective potential 

Since the T cc system is very weakly bound, a comment on long-range colour van der Waals forces 
is in order [31, 32, 33]. This forces appears due to the colour structure of the confining potential. 
They act between colour singlet clusters and their asymptotic behaviour depends on the power d 
of the long range potential. In the case of the linear confining interaction we have 

V(r) v . d . Waals = G(r d - 4 ) = 0(r- 3 ) 

This interaction appears due to the colour polarisation of two mesons in the colour singlet state. 
Such a long-range colour force is not physically allowed and is artefact of the potential approach. 
It is not present in the full QCD where quark-anriquark pair creation from the confining field 
energy would produce an exponential cut-off of this residual interaction. One might be concerned, 
however, that this spurious interaction could have some misleading effects, when the system is very 
close to the threshold, as is the case with the T cc tetraquark. To show that this is not the case, we 
present in Fig. 11 the effective potential densities 

Vij(r) = (^iVijin^Sir - nj)\il>) = V l3 {r) Pij {r). (5) 



In Fig. 11 one can see that this effect is indeed present at large separations (r > 2 fm) but it is 
extremely small. Integrating this attractive tail of the potential contributes less than 100 keV to 
the binding of the system. Another interesting feature of the effective potential shown in Fig. 11 
is the repulsion between quarks at the medium distance between quarks (1.5 fm> r > 2 fm). The 
maximal value of potential barrier is Vij(r ~ 1.5 fm) = Vij/ptj — 1 MeV, too small to produce an 
additional resonant state. 
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Figure 11: Potential densities Vij between (anti)quarks as calculated from Eq. 5 for Bhaduri 
potential. Inserted : Enlarged section of the figure, where van der Waals attraction and medium- 
range repulsion can be seen. 

At even shorter distances we have a strong attractive interaction between (anti)quark-(anti)quark 
pairs, in particular there is a strong attraction between quark and antiquark. The major part of 
this interaction bind them into the D or D* meson, since the molecular structure is dominant in 
the T cc tetraquark. The residual part of this interaction helps together with the forces between 
quark-quark and antiquark-antiquark pairs to bind the two mesons into the tetraquark. This in- 
teractions are effective at small interquark distances (r < 1 fm), where the atomic configuration 
is important. Therefore it is crucial even for the tetraquarks with molecular structure that the 
model used in the calculation is capable of describing accurately also the baryon spectra. 

D Numerics 

We solve our four body problem by diagonalization of the Hamiltonian in a space spanned by 
Gaussian function. We built our basis step by step so that at each step all configurations \ipa) 
(a = 1, ...13) are tested and parameters Cj are optimized. We then took the best configuration as 
the next base state. This procedure is very similar to the stochastic variational approach [34]. The 
main reason for using Gaussian basis is that all matrix elements can be evaluated analytically. 

We were very careful that the basis states are linearly independent so that the eigenvalues of 
the overlap matrix (ip a \ip0) 1S n °t too close to zero which would cause numerical instability. The 
dimension of the basis was between 100 (AL1 potential) and 140 (Bhaduri potential). Convergence 
of the energy of the T cc tetraquark for three different runs of the code is shown in Fig. 12. Here the 
asymptotic state of two free meson presents local minima toward which the results are converging at 
first. Only at sufficiently large number of basis states (N > 70) the bound state can be recognized. 
The initial parameters are always randomly chosen and then the optimization by Newton or simplex 
method is performed. 
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Figure 12: Energy of the T cc tetraquark with Bhaduri potential as a function of the number of 
the basis states for three different runs. The D + D* threshold is also shown. Since the initial 
parameters are chosen randomly, the convergence is similar as with the stochastic variational 
approach. 
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